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The thermal spectrum seen by accelerated observers in Minkowski space vacuum, the Unruh effect, is 
derived within the tunneling mechanism. This is a new result in this mechanism and it completes the treat- 
ment of Unruh effect via tunneling. Both Bose-Einstein and Fermi-Dirac spectrum is derived by considering 
tunneling of scalar and spin half particles respectively, across the accelerated Rindler horizon. Full solutions 
of massless Klein-Gordon and Dirac equations in the Rindler metric are employed to achieve this, instead 
of approximate solutions. 



1 Introduction 

Linearly accelerated observers should detect a ther- 
mal background with black-body spectrum in place 
where an inertial observer detects no particles. This 
result is known as the Fulling-Davies-Unruh effect 
[T] Several analysis of non-inertial observers have 
been done through a quantized field theory construc- 
tion in the coordinates adapted to an accelerated ob- 
server - Rindler coordinates. Now while analyzing 
the closely related Hawking effect f3^ for black-holes, 
an easier and more conceptually transparent single- 
particle analysis, the tunneling mechanism, was de- 
veloped in [5] . This considers a virtual pair of par- 
ticles formed just inside the horizon, one 'ingoing' 
towards the black-hole i.e. away from horizon, while 
the other 'outgoing' one traveling towards the hori- 
zon. Now, as against classical general relativity, this 
outgoing particle is taken to quantum-mechanically 
tunnel outside, with a small probability. 

The tunneling analysis has also been done in the 
Unruh effect |31 [H [71 18] where the Rindler wedges II 
and I (see Fig. [T| act as Black-hole interior and exte- 
rior regions, with the accelerated horizon playing the 
role of the black-hole horizon. However the tunneling 
mechanism could only derive the Unruh (or Hawking) 
temperature and not explicitly the spectrum. Re- 
cently this drawback was removed j9j through a quan- 
tum statistical analysis of a system of particles that 



are tunneling across the horizon, leading to a clear 
derivation of the Hawking black-body spectrum, in 
case of a spherically symmetric black-hole. Other ap- 
plications of this method may be found in [10] where 
the Kerr-Newman black hole was considerd, and in 
where black hole solutions in Lovelock gravity is 
discussed. 

In this paper, I calculate the thermal spectrum in 
Unruh effect, within the tunneling formalism. The 
methodology follows [9], but with a modification in 
the calculation of tunneling modes. Instead of the 
usually applied WKB approximations, I use full so- 
lutions of the Klein Gordon and Dirac equations in 
(3 -|- 1)-D Rindler spacetime, as the tunneling modes. 
This is possible as the flat Rindler metric is inherently 
simpler than, say, the Schwarzschild metric where full 
solution of the Klein Gordon equation is not known. 
Thus, this article provides a new and conceptually 
appealing derivation of the Unruh effect. 

2 Rindler metric and coordinate 
extension 

The path of constant, linear acceleration (say a) in 
Minkowski spacetime, is described by the hyperbola 
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Figure 1: The Rindler wedges shown on the Minkowski 
plane. pi-,T) and (x,t) are Minkowski and Rindler coor- 
dinates. The infinities lie outside the diagram, towards the 
directions shown. Tunneling occurs from wedge II to I. 

where T &i X are Minkowski coordinates, with 
X-axis being the direction of acceleration. If I con- 
sider a range (0 < < oo) of accelerated observers, 
the entire 'wedges' I and IV can be covered with the 
resultant hyperbolae. These wedges are known as 
^Rindler wedges\ A new coordinate system - the 
Rindler coordinates - can now be set up taking these 
timelike hyperbolae and corresponding spacelike lines 
(Fig. [T]) as the coordinate lines. They are related to 
the Minkowski coordinates, in wedges I and IV (indi- 
cated as subscripts), through the following relations: 



T = F(xj) sinh(ati) 
X = F{xj) cosh(afi) 
Y = yj , Z = Zj 



(2) 



where F{x)'^ = and a is a constant. The 

Minkowski Hne element ds^ = -dT^ + dX'^ + dY^ + 
dZ"^ can now be written as 



ds' 



-a^F[xfdf + F'ixfdx^ + dy^ + dz^. (3) 



This is a generalized form of the Rindler met- 
ric and for different choices of the function F{x), we 
can get all the different forms of Rindler metric seen 
in literature. Though this metric covers the entire 
Minkowski plane, the coordinates t and x change in 
the different wedges. The Rindler horizon occurs at 
F(x)^ = 0. The two wedges I and II on two sides, 
act like the black-hole exterior and interior regions 
respectively, as seen in the case of Schwarzschild so- 
lution. 

Changing to a tortoise coordinate x^, appropriate 



for the Rindler metric ([s]) through 

F'{x) 



dxi. 



^, . dx 
aF{x) 

: -lnF(x) 



(4) 



the metric ([s]) becomes 



ds' = (a 



-dt' + dxl) + dy' + dz' . (5) 



Returning to equation ([T]), though the hyperbola de- 
scribes a real accelerating particle only in wedge I, 
the hyperbolae can be analytically extended to the 
wedges II & III as X'—T' = — K. The transforma- 
tion relating the Rindler coordinate to the Minkowski 
coordinate then becomes 



T = ^(xii) cosh(aiii) 
X = F(j;ii) sinh(atii) 

y = yn , Z = Zjj. 



(6) 



Now the maximal extension of the Rindler coordinate 
system (t, x, . . .) is the Minkowski (T, X, . . .) which 
is defined everywhere throughout the four Rindler 
wedges, i.e. the entire Minkowski plane, irrespec- 
tive of accelerated horizons which the Rindler ob- 
server encounters. The Rindler coordinates however 
undergo a finite shift through the horizon, as can eas- 
ily be seen through a comparison of ^ and ([6]). A 
relation between the t — x^, coordinate pair outside 
(Region I) and inside (Region II) can then be writ- 
ten as 



til — 



ITT 

2a 
in 

2a' 
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The y and z coordinates on the other hand, remain 
unchanged across the horizon. It is to be noted that 
the transformation = + i^! = ~ 1^ 

suffice in relating the coordinates II and I. However 
this second pair leads to some problems in taking the 
classical limit of the tunneling probability, and so is 
not considered at the accelerated horizon, as will be 
discussed later. Such transforms were reported ear- 
lier in da [13]. 

3 Wave function: Scalar particles 

The massless Klein Gordon equation = 0, writ- 
ten in the Rindler metric ([5]), reads 

p—2aXi, 

{-d'^^ + dl^^-\ + + dl^ = 0. (8) 
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Since the metric is independent of the coordinates 
t,y Sz z,l take an ansatz for $ as 



(9) 



where 17 is a constant. This O is related to the 
locahy observed energy iv at some Rindler space- 
time point through a red-shift |14| [T5] relation 
El Vi = E2 V2 = ^ connecting the observed ener- 
gies El & E2 at two different points in a gravitat- 
ing system at equilibrium. The result ensures that 
though the observed energies E and the Tolman red- 
shift factor V vary locally as functions of the coordi- 
nates (here x^), their product 17 is a constant. In a 
gravitational system in equilibrium, this condition of 
the constancy of J7 characterizes the equilibrium, just 
as is done by temperature in a laboratory thermo- 
dynamic system at thermal equilibrium Here, 
this locally observed energy E is the energy uj of 
the tunneling-particle and the quantity V becomes 
Y^l^ool- Thus I have 



ijj ae 



a 



(10) 



where appropriate definitions -F(x^) = 1/a = e"^* 
(see metrics [3] & [5]) have been used. Substituting the 
ansatz ([9| in the Klein Gordon equation ([s]), yields 
the following differential equation in 0(x^) 



A-2 



0(Xv, 



0, 



(11) 



U2 I U2 



with = 

Some observations can immediately be made from 
equation (11). Near the horizon, as x* — > —00, the 



term containing k± drops out and a simple harmonic 
type equation with plane wave solutions is obtained. 
Again at large spatial distances, x* — > 00, and now 
the term containing $7^ becomes negligible. This 
leaves an equation with exponentially increasing and 



h ^ 



and K{ 



decreasing solutions /q 
where /q and Kq represent the zero-th order mod- 
ified Bessel functions of first and second types re- 
spectively. Thus throwing away the Iq solutions, we 
have an exponentially vanishing solution at infinity, 
in Kq. Similar conclusions have also been reached at 
by Boulware fTB]. 



The solution of the full equation (11) that is well 
defined through the horizon is. 



4){x^) =A. 



TVSI 

e 2afi T 



" 2ah T 



1 



1 + —] Lri 
an J ah 



+ 



(12) 



where Azp are arbitrary integration constants. For 
small arguments, the appropriate expansion of the 
modified Bessel function is Iu{z) ~ p^j^^. This 
holds if /c_L <C 17 and also especially near the hori- 
zon where x^, - 



00. Therefore (12) simplifies to 



2h 



+ aft 



The total wave function <I>(t, x^, y, z) near the horizon 
is then 

<^{t,x^,y,z) =B.^e-Tii^^^+^*^+''yy+''^'^ + 

Soute"s[f^{*-^*)+fcyf+'=-^l, (13) 

with all the constants clubbed together within 
n/out- The subscript 'in' here stands for the ingoing 
mode which travels toward the accelerated horizon 
at X* = —00, while the subscript 'out' stands for the 
outgoing mode traveling away from horizon, i.e. to- 
wards Xj, = 00. 

4 Wave function: Spin | particles 

Spinors are introduced on a general curved space- 
time with metric g, by going to a local Lorentz frame 
with metric r], at each spacetime point |T7]. This 
procedure is best done by constructing tetrad fields 
which map curved spacetime tensors to local Lorentz 
frame and vice-versa. For a Rindler metric in the 
tortoise-like coordinate system the tetrad field 
is defined through the relation 

= y' vab 

where latin (a, 6, . . .) and greek (/i, z^, . . .) letters run 
over local Lorentz and curved space indices respec- 
tively. The explicit choice of the tetrad field 
adopted here is 



= diag(ae''^%ae''^%l,l) 



(14) 



and the metric signature, both global and local, is 
kept same as (— , -|-, -|-, -|-). 

The massless Dirac equation is written as |18| |2] 



(15) 



where 7*^ are the Dirac matrices obeying the usual 
algebra [7^,7^] = ^rj""^ and F^ are connection coeffi- 
cients given by 



F 



1 



ab 



2 

1 - 

4 _' 



bu]ij, 
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The covariant derivative over the curved space in- 
dex of the tetrad is defined in the usual way Vbu-,^ = 
d^Vhu — T'^i^Vba, where F^^^ is the Christoffel symbol. 
On using the properties of 7 matrices and the diag- 



onal choice of the tetrad (14), the spin-connection 
becomes T, 



2 " " 

tion (|15|) then turns out to be 



1 S'*'' r^, Vbx. The Dirac equa- 



—ae 



0. (16) 



The ansatz for the spinor ^ is taken as 



■0(2^* 








(17) 



Upon using this ansatz, equation (16) can be cast 
into a Schrodinger like equation 



where the Hamiltonian-like operator H is 



(18) 



(fc,7V + A:.7V). (19) 



The above equation, on squaring, gives H'^il) = 0^^/^. 



Now using the ansatz (17) and adopting the con- 
vention for the gamma matrices as 7*^ = ( "^^ ^ ) , 

7-' = cT'' ) (^^^^^ j ~ 1;2,3 and are the 

Pauli matrices), the following equation for the spinor 
component functions is obtained 



n? 



♦ = 0. (20) 



Here ♦ stands for the functions ^(x^,) and B{x^). As 



in the case of scalar particles (11), a study of asymp- 



totic behaviour of this equation show that solutions 
near the horizon are oscillatory, and that near to in- 
finity are vanishing in nature. Solution for the full 



equation (20) turns out to be 



m1*^ r ( 1 
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(21) 



It is to be noted that both (21) and (12) are full 



solutions of the respective differential equations ( 20 ) 



and (11 ), without any approximations of parameters. 



This can be easily verified by substituting back these 
solutions in the original differential equations, to see 
that they satisfy them without using any approxi- 
mations on the parameters. However, an approxi- 



mation is used to get a simpler solution from (21), 



by using the appropriate expansion of the modified 
Bessel function for small arguments which is given 
as Iu{z) ~ rli+u) • This holds if A;_|_ ^ and also 
especially near the horizon where x^ —00. So, for 
the region near to the horizon, the total spinor ^ can 
finally be written as 



^ g-j[[^{t-x*)+kyy+k^z] /22^ 



with i^in/out being constant spinors and the subscript 
'in' or 'out' standing for 'ingoing' or 'outgoing' modes. 



Unruh effect through tunneUng: 
thermal spectrum and tempera- 
ture 



Uptil now, I had found single-particle wave-functions 
for bosons and fermions, by solving for the Klein- 
Gordon and Dirac equations in the Rindler coor- 
dinates. These solutions are valid in both Rindler 
wedges I and II, but in coordinates {tj_,Xi) and 
(tiijXji) respectively. Classically, both the ingoing 
and outgoing modes (say of a virtual pair instanta- 
neously produced) in wedge II are trapped, as noth- 
ing from inside can cross the horizon and come out 
to wedge I. However in the tunneling mechanism, an 
outgoing particle can quantum-mechanically tunnel 
out across the horizon and into wedge I. This process 
occurs with a probability given by the Maxwell term 
e fta , that appropriately goes to zero in the classical 
(^ — > 0) limit. Now to find the energy distribution 
of a collection of such particles, I will (following [5]) 
construct a suitable density matrix for both bosons 
and fermions, and find out the average number of 
particles having some particular energy uj. 

Starting first with bosonic particles, the relation 
between inside and outside wave-functions is found 
by using the connection between coordinates ([T]) in 



4 



equation (13) for the modes. 



^-j^[fl(tii+Xi,jj)+kyyii+k^zjj] _ 

g-^[n(ti+a::*i)+fcyj/i+fcz2i] 
p „-T[f2(tii-a;*ii)+fcayii+A:z2ii] _ 

(23) 

Now, let there be 'n' pair of free particles (ingoing 
and outgoing) in wedge II. The total state of this 
system of particles, with each being described by the 



sector II modes in equation (13), is 



n=0 

oo 



{e~^\ K) |n° 



(24) 



n=0 



where Nb is a normalization constant defined 
through {xb\Xb) = 1- The sum over n runs from 

to oo here, in the case of bosons. But in case of 
fermions, as will be used later, n is limited to and 

1 by Pauli's exclusion principle. The normalization 
of \xb) leads to 

oo 

o ^ — > (n+m)7rn / ■ \ / ■ \ 

Nl e — [{mT\ {m["\j (^K) K"')) = 1 

n,m=0 

-1 



E 

.n=0 



(25) 



and finally for bosons, we have 
Nb 



2Tvn \ 2 

1 — e ah 



(26) 



The density matrix operator for this system of 
bosons is defined as usual 



Pb 



_ 2ttS2 
1 — e a,h 



E 

n,m=0 



(n + m)7rf2 

e ^ |n") \n°" 



(27) 



Since ingoing waves are trapped within the horizon 
and outgoing particles contribute to spectrum, we 
trace out over the ingoing particles, to form the den- 
sity matrix for outgoing modes. 



>T = - e-^]Y,e- — \nr){nr\. (28) 



n=0 



The average number of outgoing particles is then cal- 
culated as 



n=0 

1 



27rn 
e ah - 

1 



1 



2TTU1 

e 



1 



(29) 



where in the last step, the red-shift definition of 



given in (10) was used. This is immediately recog- 



nizable as the Bose-Einstein distribution for a black 
body at a temperature T^, the Unruh temperature 
[I], given as 



ha 
2ti 



(30) 



with a being the local acceleration. 

For fermions, the same method goes through step 
by step. The connection between spinorial wave- 
functions in wedges II and I is obtained by using 
(0 in (|22). 



^ g-|;[0(tlI+Xtii) + fcy?/ii + fcz2ii] _ 

^ g-|;[r2(tl+a;*i) + feyyi + fcz2i] 

t -j[0.{tjj-Xi,ij)+kyyjj+kzZij\ _ 

Sout^ 

(e" ^ ) Lu.e- s m^-^*^)+kyyi+k.zi] 

(31) 

The normalization of the total state ket for fermions 

1 



\Xf) 



n=0 



n=0 



(32) 



is again done through {xf\Xf) = 1- The sum over 
number of particles in a given state, in fermionic cal- 
culations, always run from to 1, following Pauli's 
exclusion principle. The normalization constant Np 
turns out to be 



_ 27rn 
1 + e ah 



(33) 



The density operator for fermions pp is defined as 
\Xf){Xf\- Using the outgoing fermionic density oper- 
ator pPp^ the spectrum is calculated as average number 



5 



of outgoing particles 



{hp) 



1 

e ah + 1 

1 

e fia + 1 



1 + e aft „=o 



(34) 



where equation ( 10 ) was used in the last step. This is 



the Fermi-Dirac distribution, at the Unruh temper- 



ature Tjj defined in (30), and constitutes the Unruh 



effect for accelerated fermions [19j. 

6 Discussions 

In this paper, I calculated the thermal spectrum for 
Unruh effect within the tunneling mechanism. The 
Unruh temperature was identified via a comparison 
between the calculated thermal distribution and the 
standard form of Bose or Fermi distributions. 

However, the temperature can also be calculated 
directly via tunneling. Say, when an observer in 
Rindler wedge I observes an outgoing particle, com- 
ing from within wedge II, he will see the wave func- 
tion from inside wedge II with a factor as shown 



in equation (23) for scalar particles, and in (31) 



for fermions. The ingoing wave, however, does not 
change by such a factor between the two wedges. 
To calculate the temperature, we can now use the 
principle of detailed balance, = e~T, where 

^in 



P. 



out /in 



I wa ve- function p is the outgoing/ingoing 
probability, w is the observed energy at that point 
and T is the temperature. This gives the Unruh tem 
perature Tu 



^ after suitably using the red-shift 



definition of given in (10) at the observer's point. 



For a more detailed discussion on the connection be- 
tween the earlier approaches to tunneling and the one 
used in this paper, the reader is directed to [9, ilS] . 

Another point to note concerns an apparent am- 
biguity in the sign of the factors chosen to con- 
nect the coordinates between wedges I and II in 
equation ([7]). It can be verified that the relations 
tji = ti + ^'i = x-ki — 1^, which have a change in 
sign between the t and term signs, also connect the 
transformations ^ and ([6]). However this set is not 
employed as it produces an exponential factor with a 

wrong exponent sign { eaU ] in the relations between 



wave-functions in wedges I and II (23 & 31). With 



this alternate sign, the probability of outgoing parti- 
cles to tunnel out across the horizon, diverges at the 



classical limit of ^ ^ 0. This is clearly unacceptable. 
So this obvious criterion of conformity with classical 
limit is used to remove the above said ambiguity of 
signs. 
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